Mon. Not. R. Astron. Soc. 000,[T]l7](2009) Printed 3 September 2009 (MN KTeX style file v2.2) 



A Mechanism for the Oxygen and Iron Bimodal Radial Distribution 
Formation in the Disc of our Galaxy 
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ABSTRACT 

Recently it has been proposed that there are two types of SN la progenitors - short-Hved 
and long-lived. On the basis of this idea, we develope a theory of a unified mechanism for 
the formation of the bimodal radial distribution of iron and oxygen in the Galactic disc. The 
underlying cause for the formation of the fine structure of the radial abundance pattern is 
the influence of spiral arms, specifically, the combined effect of the corotation resonance and 
turbulent diffusion. From our modelling we conclude that to explain the bimodal radial dis- 
tributions simultaneously for oxygen and iron and to obtain approximately equal total iron 
output from different types of supernovae, the mean ejected iron mass per supernova event 
should be the same as quoted in literature if maximum mass of stars, that eject heavy ele- 
ments, is 50 A/0. For the upper mass limit of 70 A/0 the production of iron by a supernova II 
explosion should be increased by about 1.5 times. 

Key words: Galaxy: abundances - Galaxy: structure. 



1 INTRODUCTION 

' Galactic nucleosynthesis studies provide a very important informa- 
tion about the structure and evolution of the Milky Way Galaxy, 

. since chemical abundances accumulate and retain signatures of 
milestones in its history. This is why this research area is often 
called the "cosmic archaeology". Knowledge of heavy element en- 
richment in galaxies plays a crucial role in various scientific areas. 

One of the most prominent features of the chemical compo- 
sition both in the Milky Way and other galaxies is the metallicity 
gradient. For a long time, it was widely believed that the radial 
metallicity distribution in galactic discs is described simply by a 
linear function, with the more or less consta nt gradient value be ing 
typical for the most part of the disc (see e.g. lwielen et al.ll 19961). 

Such a n oversimplifi e d bel ief was first broken by Shaver et al.l 

( Il983h and iTwarog et al.l jl997t) . They concluded that the abun- 
dance distribution in the disc of the Galaxy is not described by 
a linear function with a constant slope . In a seri e s of papers, 
lAndrievskv et alj ^20021 hereafter A02). iLuck et alj j2006l here- 
after LKA06) and others, have shown that in the Milky Way, oxy- 
gen and iron demonstrate a bimodal distribution along the galactic 
radius. In the inner part, from about 4 to 6.6 kpc, the gradient is 
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quite steep (for iron it is of the order of -0. 13dex/kpc), whereas 
in the outer part, approximately from 6.6 to 10.6 kpc, the distribu- 
tion is plateau-like (the gradient is about -0.03 dex/kpc, which is 
about 4 times smaller in absolute value than in the inner part), the 
distributions being similar both for oxygen and iron. 

This result has posed a very difficult problem that remained 
unsolved for several years. Indeed, at the epoch, when most re- 
searches believed that the abundance radial distribution is linear, 
they did not pay much attention to the fact that the gradients show 
up in observations of different elements. They also did not worry 
about that oxygen and iron have close radial distribution in the disc 
of the Galaxy. The reason perhaps was as follows. The radial distri- 
butions of both types of SNe, averaged over the azimuthal angle in 
the galactic plane, are approximately exponential with close radial 
scales. Hence, we may expect that the abundances, expressed in the 
logarithmic scale, will be approximately linearly distributed along 
the galactic radius with close gradient values. But the bimodal ra- 
dial abundance pattern breaks this oversimplified representation. 

In the present p aper, we continue the res earch started by 
iMishurov et al.1 ( l2002h and lAcharova et al.l feOOSi) . In the cited pa- 
pers and others a theory of formation the bimodal abundance pat- 
tern under the influence of spiral arms was developed. The main 
cause for formation of the bimodal distribution is assumed to be the 
corotation resonanc e which is located close to the solar galacto- 
centric distance (e.g. iLepine et al]|200lL and papers cited therein). 
In the previous papers we have only considered the enrichment of 
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the interstellar medium (ISM) by oxygen, because the main source 
of oxygen - SN II - are strongly concentrated in spiral arms. Their 
progenitors are massive stars with extremely short life-time, thus, 
SN II do not escape from spiral arms where they were born. Hence, 
oxygen is the most obvious indicator for the arm influence on nu- 
cleosynthesis (see details in the papers cited above). Our earlier 
results have strongly suggested that other models of spiral arms, 
like recurrent wave patterns with differ ent locations o f the corota- 
tion resonance during galactic life (e.g. ISellwood & Binnev 20o3) 
or wi th the corotation in the inner part of the Galaxy I Weinbera 
1 1994 . etc., failed to explain the bimodal radial distribution of oxy- 
gen in the galactic disc. 

Some other explanations have been proposed recently for the 
decrease in abundances of many elements with galactocentric dis- 
tance as well as for the bimodal structure of the observed gradient. 
These explanations are mostly centered on the radially depen dent 
infall rate and/or star formation efficiency (see, e.g., Colavitti et al.l 
I2OO9I ; IFu et ai]|2009l ; iMagrini et alj|2009l) . In various models the 
gradient may have both internal (inside-out disc formation) and 
external (dilution by metal-poor infall from the halo) origin with 
respect to the disc. In particular, IColavitti et al.l ( 1200 9*) concluded 
that the inside-out disc formation and the density threshold are both 
necessary ingredients for a model which is supposed to reproduce 
simultaneously all the disc constraints, including steep gradient in 
the inner disc and flat gradient in the outer disc. However, in studies 
cited above the gradient slope changes at 11-12 kpc, while obser- 
vations indicate that the change occurs closer to the Galactic centre, 
at around 7 kpc (LKA06). 

As it was noticed above, the observed radial distribution of 
iron is similar to oxygen. Unlike oxygen, only ~ 1/3 of iron is 
produced b y SN II. About 2/3 of it is synthesized during SN la 
explosions ( lMatteuccill2004 . For a long time, precursors of these 
events were considered to be quite old objects with ages of the order 
of several (4—5) billion years. Even if they are born in spiral arms, 
for a long time prior to the explosion they "forget" the initial affinity 
to spiral arms due to both azimuthal and radial migration. So, at first 
sight, we should not expect strong manifestations of galactic spiral 
arms influence on radial distribution of iron. Hence, the similarity 
of bimodal radial distributions of oxygen and iron in the disc of the 
Galaxy repres ents a difficult enigm a. 

Recentlv. lManucci et alj l [20()6h have shown that there may ex- 
ist two sub-populations of SNIa. Events of the first type (short- 
lived progenitors, they were called "prompt") have ages less than 
10* years, while events of the second type ("tardy") have long- 
lived progenitors, with delays between formation and explosion be- 
ing measured in bi llions years (the theor etical basis for such kind 
of division see in ilben & Tutukovl jl984 ). The frequencies of ex- 
plosions and iron outputs for both populations are approximately 
equal. Prompt SNe la precursors, being quite massive stars, retain 
their relation to the spiral pattern. Taking into account that SN II 
and short-lived SN la jointly produce ~ 2/3 of iron, this new result 
opens an opportunity for solving the above problem in the frame- 
work of the theory of spiral arms influence on radial abundance 
distribution. 

In this paper, we propose a unified mechanism for formation 
of bimodal radial distribution both for oxygen and iron. 



2 OBSERVATIONAL DATA 

In Fig. □ the radial distribution of {[O / H]) and {[Fe/H]) 
vs galactocentric distance r is shown (as usual [X/H\ — 
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Figure 1. Radial distributions of oxygen (upper panel) and iron (low panel) 
averaged over the azimut hal angle in the ga l actic p lane. Open squares are 
the observed values from lAndrievskv et alj j2002l) and iLuck et al] j20O6h 
averaged within bins of 0.5 kpc width. The error-like bars are the stan- 
dard deviations of the corresponding mean values. The abundances of A02 
and LKA06 were co rrected for the new solar composition according to 
lAsplund etai] )2006l) . Position of the Sun on the diagrams is shown by the 
usual solar mark. 



log{Nx/NH)star ^ log(iVx /A^H ) , where Nx,h is the num- 
ber of atoms for X element or for hydrogen, correspondingly). 
We divided r into bins of 0.5 kpc and averaged the correspond- 
ing data from A02 and LKA06 within each bin. The angle brack- 
ets ({. . .)) denote these averaged values. The error bars show the 
standard deviation of the averaged values within the bin. We used 
the solar galactocentric distance = 7.9 kpc, as in LKA06, 
but abundances wer e corrected for the new solar composition 
jAsplund et"ai]|2006l) . In this figure, the Sun is shown by the usual 
mark. The similar bimodal patterns in oxygen and iron radial dis- 
tribution are clearly seen in Fig.[T] 

Our aim is to explain the bimodal radial distribution of oxygen 
and iron in the most part of the Galaxy, i.e. the steep gradient in 
the inner part of the Galactic disc (r ^ 7 kpc) and rather shallow 
distribution up to about 10 kpc with a bending at approximately 7 
kpc. The peculiarities in the distributions at r > 10 kpc are beyond 
the scope of our theory. 
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3 BASIC EQUATIONS AND MODEL 

T he equations of chemical evol u tion of ISM are clos e to the ones 
of lTinslevI JTgSOh . lChiosil ([T980h . lLacev & Falll (Il985h : 



P-g = f - i' + / (m - m^)ip(t - T^)(t>{m) dm, 



(1) 



(m — rriji,) Zi(t - 



^)(t){m) dm + 



r or \ 



dr 



(2) 



where jig^i are surface mass densities for interstellar gas (subscript 
"g") and heavy elements (where subscript "i" is "O" for oxygen 
and "Fe" for iron), Zi = /x^/ /ig is the mass fraction for the cor- 
respondin g elements, tj} is the star formation rate (SFR), for which 
we adopt (lKennicuttl[T998l) : 



(3) 



where v k & normalizing coefficient. All the above quantities are 
functions of time t and radius r, but we do not explicitly indicate 
them except for the cases when there is the time shift, t — r^, in the 
integrals, m is the stellar mass (all masses are in solar units), jtil 
and mu are the minimum and maximum ste llar masses, =0.1 , 
while for mu we use two values 50 and 70 jTsuiimoto et al.ll995h . 
Our calculations show that the results are very close for the both 
above values of mu- That is why in our figures we only show the 
radial distributions for mu = 70. The life-time of a star of mass 
m, Tm (in Gyrs), on the ma in sequence is given by \o^(Tm) = 
0.9 - 3.8 log(m) + log^(m) ('Tutu kov & Krugelll980l) . The initial 
mass function <j>{m) is Salpeter's one with the exponent of - 2.35. 

In full formulation of the problem, we should take into ac- 
count the divergent terms, which describe the influence of radial gas 
flow ("radial inflow") wit hin the ggJactic disc. It s effect was consid- 
ered in various models by'Lacev & Fall (19851. |Portinari & Chiosil 
( |2000|) .I Mishurov et al .l 0002). In the present paper, we will not in- 
clude the radial redistribution of the interstellar gas. This problem 
will be separately considered in details elsewhere. 

The infall rate, /, of intergalactic gas onto the galactic disc we 
describe as: 



/ = B(r)exp(- 



t_. 



For the time scale t / the following representations were used: 
l)t f = 3 Gyrs (let us call this law as "rapid" disk formation model, 
see iPortinari & Chiosil h999l) ): 2) tf = 7 Gyrs ("slow" disc for- 
mation) and 3) tf — c + dr, where c and d are some constants. 
The last dependence produces the so-called "inside-out" scenario 
for the galactic disk formation. W e consider various va lues for c 
and d. i) c = 1.033; d = -1.267 jChiappini et"Zll200lh . It is ob- 
vious that this representation may not be used at r = 1.226 kpc 
since here / ^ 0. To exclude this peculiarity we use the above 
representation forr>r/ = 2 — 3 kpc. Inside r/ let us suppose 
tf = a + brf = const. Our experiments show that in the range of 
galactocentric distances, interesting for us, the res ult does not de - 
pend on the exact value of r/. ii) tf = 1 + 0.76r i IFu et al.ll2009h . 
At the Sun location the both expressions give tf = 7 Gyrs. 

To fix B{r) as usually we adopt that the final radial to- 
tal density distribution (stars -I- gas) at current epoch, t = Td 
(Td = 10 Gyrs is the disc age) is exponential with the radial scale 



rd = 2.5 - 4.5 kpc (see, e . g..lNaab & Ostrikeil ^20061) , IFu etal] 
( l2009l) , ISch6nrich & Biiiiievl J2OO8I) and others). 

In all our experiments we normalize the infall rate so that the 
total present-day density at the solar location to be equal to 50 Mq 
jNaab & Ostrikej|2006l) . To obtain the normalizing coefficient v 
for SFR function (see Equation (3)) we suppose that at the solar 
galactocentric distance /Xg « 10 Mq. 

The infalling gas is assumed to have abundances Zi^f — 
0.02 Zi,0. 

The value for the mass of stellar remnants (white dwarfs, neu- 
tron stars and black holes) is adopted as follows: for m ^ 10 
m» = 0.65m'''^^^; in the range 10 < m < 30 m» = 1.4; if 
30 ^ m < mu the remnant is a black hole with m„ = 10; fi- 
nally for m ^ mu the stars are black holes right away from their 
birth and they are r emoved from the nucleosynthesis and return - 
ing the mass to ISM jMassevich & Tutukovl ( ll988l) . lBretor] ( |2009|) . 
IWooslevetai] ( l2002h ). 

The last term in Equation (2) describes the radial diffusion 
of heavy elements. We need it to smooth out strong depression 
in abundance distribution forming when the role of spiral arms on 
heavy elements synthesis is taken into account (in the case of low 
gradient, the role of the diffusion happens to be negligible). Note as 
well, that the diffusion smoothes out the elemen t distribution as ad- 
mixtur es but not the gaseous density as a whole dLandau & Lifshit3 
l[l98i)). The nature of the diffusion we attribute to turbulence in 
ISM arising due to stellar winds, supemovae explosions, galactic 
fountains, etc. For a phenomenological description of this process, 
we introduce the diffusion term into Equation (2). To estimate the 
diffusion coefficient D we model the turbulent ISM by a system 
of clouds, so a simple gas kinetic approach enables us to derive th e 
expression for D (see iMishurov et al .l2002l ; lAcharova et al.l2005h : 



D 



Hvt 



where H = 130 pc is the thickness of the gaseous disc, vt ~ 
6.6 km s"'^ is the one dimensional dispersion of cloud chaotic 
velocity, m^ and a are the mass and rad ius of a typical cloud, 
mc/a^ = 100 Mq pc"^ ( lElmegreenlll987h . 

Enrichment rates of ISM by SNIa and SNII, correspond- 
ingly, are set by values Ej'^ and Ej^. The rates of oxygen and 
iron synthesis in SN II events are described by similar expressions 
Ej^ — r]Pj^R^^, where P/^ is the mass of ejected oxygen or iron 
per one SNII explosion, R^^ is the rate of SNell events 



R'\r,t) 



(1 - A) / t 



Tl)(t){m) dm. 



ms — 8, A — 0.0025. The last parameter gives the fraction of 
binary systems pr oducing SNe la relat ive to the whole mass range 
jManucci et aLll2006i) . The factor r] describes the 
influence of spiral arms (see below). 

Contribution from SN la to the enrichment is represented 
as EI^ = Ei^p + Ei^T since there are two populations of 
SNe la progenitors, namely, short-lived ("P ") and long-lived ("T", 
iManucci et al.ll2006l ; lMatteucci et al.ll2006h . As we wrote in Intro- 
duction, being massive stars, the short-lived progenitors are con- 
centrated near spiral arms. So, by analogy with SN II rates we write 
down Ei,p = rij Pl^E^p , where 7 is a correction factor. For long- 
lived SN la Ei^T = C P/^-Rt • Here P^^y are rates of SN la ex- 
plosions for P- and T-progenitors, P}^ is the mass of ejected i-th 
element per one SN la event, and C, is another normalizing coeffi- 
cient. 
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To compute -Rp y we follow iMatteucci et alj ( |2006|) but ex- 
plicitly separate the "Delay Time Distribution" function, DTD{t), 
into two parts: DTD = Dp for t ^ ts = 10"^ °'' Gyrs (P- 
progenitors): 

log(Dp) = 1.4 - 50[log(r) + 1.3]^ 
and DTD = Dt for r > rg (T-progenitors): 

log(DT) = -0.8 - 0.9[log(r) + 0.3]^ 

So, 

R^pir, t) = Akc, j ^{r, t - T)Dp{T)dT, 

Tg is the life-time for a star of mass m = 8 and 

t 

Rriryt) = ylfcc / ^/.(r, t - T)DT{r)dT, 



200 



ka = J 4>{^) dm (see details in lMatteucci et al.ll2006h . 

mi 

Notice that Dp has a very sharp maximum at r = rp = 
10-^-^ Gyrs « 50Myrs. This means that for rg ^ t < rp the 
main contribution to i?p comes from the vicinity of r ~ t whereas 
for t^Tp the vicinity of rp is most important. The corresponding 
integral was computed asymptotically using the Laplace method. 

To incorporate the influence of spiral arms on formation of ra- 
dial abundance pattern in our theory, we follow the idea o f' Oorll 
( Il97 4). For this, we take into account that SNII and short-lived 
SNIa progenitors are concentrated near spiral arms. Hence, the 
enrichment of any volume of ISM by heavy elements only takes 
place when the volume is close to the nucleosynthesis sites, i. e., 
inside or near spiral arms. The more frequently the volume passes 
an arm, the higher the rate of its enrichment gets. The frequency 
of the volume entering a spiral arm is proportional to the differ- 
ence \^{r) — Q.p\ where Sl(r) is the angular rotation velocity of 
the galactic disc, f2p is the rotation velocity of the wave pattern 
responsible for the spiral arms. So, we write r; = P\Q.{r) — Q,p\Q, 
where /3 is a normalizing coefficient, G is a cutoff factor. The cor- 
rection factor 7 in expression for Ei^p was introduced to take into 
account that we aim to obtain close amounts of iron produced by 
SNII and prompt SNIa. 

We assume that long-lived SN la progenitors do not concen- 
trate in spiral arms and are uniformly distributed over the azimuth 
angle in the galactic disc. Therefore, unlike r;, C = const. 

Spiral arms exist between inner, r^, and outer. To, Lindblad 
resonances, their locations being given by the condition Q.{ri) — 
K,{ri)jn ^ Sip ^ Q.{ro) + K,{ro)/n. Here, n is the number of spi- 
ral arms, k is the epicyclic frequency. So in the wave zone, i. e., be- 
tween Lindblad resonances, G = 1, and beyond th em G = 0. From 
the galactic density wave theory jLin et al.llT96^ it is known that 
the density wave pattern rotates as a solid body, i. e..,Q,p = const 
whereas the galactic matter rotates differentially: is a function 
of the galactocentric distance r. The radius rc where both veloci- 
ties coincide {Q.{rc) = Sip) is called the corotation radius. In the 
absence of diffusion from the above representation for t] one can 
expect the formation of a valley (or a gap) in the abundances distri- 
bution near the corotation (see below). The diffusion smoothes out 
the valley so that the combine effect of it and the corotation reso- 
nance results in formation of the bimodal radial distribution of the 
abundances. 
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Figure 2. The rotation curve Sl(r) and its combina t ions w ith epicyclic 
frequency k approximated over the data of lClemenj il985() adjusted for 
tq = 7.9 kpc. The horizontal dashed line shows the adopted value of 
the angular rotation velocity, Hp , for the density wave pattern which cor- 
responds to the location of the corotation resonance at ~ 7 kpc. The 
vertical dotted lines indicate the positions of Lindblad resonances. 



In otir mod elling we used a rotation curve based on CO data 
dClemensll 19851) but adjusted for tq = 7.9 kpc: 



,3.6, 



0.1, 



rSl = 260 • exph— - { — )'] + 360 ■ exp(-— - — ). (4) 
150 r 3.3 r 



The rotation curve is shown in Fig.|2] Also indicated in Fig. [2] are 
locations of Lindblad and corotation resonances derived for Sip = 
33.3 kms^^ kpc^^ (correspondingly rc ~ 7 kpc). 

Mathematically our system of equations splits into two 
groups. Equation (1) is an ordinary differential equation. Solving 
it, we find both the gaseous and stellar density radial profiles at any 
time t, fitting the observed stellar and gaseous densities to the ones 
at present time. Equations (2) are differential equations in partial 
derivatives. For the boundary conditions we adopt the natural ones 
at the galactic centre and at the galactic peri phery which guarantee s 
the finiteness of our solutions (see details in lAcharova et alj2005l) . 

Corresponding constants for "chemical reactions " (/3, 7, Q 
were derived in the following way. As can be seen from lMatteuccil 
( l2004l) . lTsuiimoto et aO ( Il995h SNe la produce a small part of oxy- 
gen, no more than several percents. Hence, at first step we can ne- 
glect by oxygen enrichment from SNe la. This enables us to find 
independently the parameter (3 by fitting the observed radial distri- 
bution of oxygen. After that we derive other two quantities, 7 and 
C,. To do that we impose the additional condition that short-lived 
(P) SN la, long-lived (T) SN la, and SN II produce equal amounts 
of iron (~ 1 /3 each). This constraint enables to solve the problem 
unambiguously. The above derived values for chemical reactions 
constants we consider as zero order approximation for them. After 
that, we include the contribution to oxygen production from SN la 
and find the revised chemical constants. 
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Figure 3. Theoretical abundance patterns derived for constant tf =3 
Gyrs (rapid disc formation model) superimposed on the observations. Dot- 
ted lines are for = 2.5 kpc. Solid lines - for =3.5 kpc. Dashed lines 
- for ra = 4.5 kpc. 

4 RESULTS AND DISCUSSION 

Below we show the results of our modeling for the adopted value of 
Q p which corresponds to the location of the corotation resonance 
atrc ^ 7kpc.E] 

In Fig. [3] and Fig. |4] are shown the computed abundance pat- 
terns for constant tf. In the case of rapid disc formation (tf — 3 
Gyrs) we demonstrate the dependence of the distributions on the 
radial scale r^. It is seen that in the inner part of the disc the corre- 
sponding curves are indistinguishable. They slightly differ beyond 
the solar location. Overall, the coincidence of the model with the 
observations is good but perhaps it is the best for the radial scale 
ra — 3.5 kpc. 

In the case of slow disc formation (tf = 7 Gyrs, Fig.|4l( we 
failed to fit the model to the observations. The matter is, by means 
of one set of chemical constants we can fit the observations, say, in 
the inner part of the disc, but for matching with the data in the outer 
part we have to use another set of the constants. 

The general conclusion is: in the framework of scenario of 
disc formation with constant tf we can succeed in agreement with 

^ The above value for the corotation radius is not to be considered as a 
determination of Tc, it only produces a good lit to the bimodal radial dis- 
tribution of oxygen and iron in a simple model that does not include all the 
fine effects. 
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Figure 4. The same as in Fig. 3 but for tf = 7 Gyrs (slow disc forma- 
tion model) and = 3.5 kpc. Dotted lines are for the set of the fitting 
parameters fi = 0.030, 7 = 0.32, C, = 0.22. Dashed lines are for the set 
P = 0.037, 7 = 0.30, C = 0.25. From this figure is seen that there is 
no way to reconcile the model with the observational data in all range of 
galactic radius by means of the same set of chemical constants. 



the observed radial distributions both for oxygen and iron in the 
galactic disc with short time scale ~ 3 Gyrs. The fitted param- 
eters happen to be as follows: for the upper stellar mass mu = 70 
(3 = 0.019 , 7 = 0.44, C = 0.18; for m[/ = 50 /3 = 0.026, 
7 = 0.32, parameter C, remains the same. 

However, to derive the observed distributions and satisfy the 
condition that the above three types of sources produce approxi- 
mately equal amount of iron, we had to adopt the following masses 
ejected per one SN. In the c ase mu = 50 the corre sponding val- 
ues are equal to the ones in iTsuiimoto et al.] l ll995h , i. e., Pq = 
1.80 Mq, Pll = 0.613 M0, Pll = 0.084 M0. For the upper 
value of mu = 70 Pq = 2.47 Mq, P^l = 0.613 M©, like in the 
last cited paper but PpJ = 0.13 M©, (i.e. about 1.5 times greater 
than in lTsuiimoto et al.lll99^ . 

Our experiments show that SNIa produce about 2% of oxygen. 

In Fig. [5] are shown the computed distributions for the 
inside-out scen ario (since they are close both for the models of 
IChiappini etal] (^2001) and I Fu et al.l f2009.) we only show the re- 
sults for Chiappini-like representation of tf). In our opinion the 
coincidence of the theory with the observations both for oxygen 
and iron is excellent. The chemical reactions are close to the ones 
of constant t/, namely: for mu = 70 /3 = 0.021 , 7 = 0.46, 
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Figure 5. The same as in Fig. 3 but for inside-out disc formation scenario 
for Chiappini - like model (see Section 3) and = 3 kpc (in the radius 
range interesting for us the results slightly depend on the exact value of rj). 
Other designations are the same as in Fig. 3. 

C = 0.18; for mc/ = 50 /3 = 0.032, 7 = 0.30, parameter C, again 
remains the same. 

To feel the role of the turbulent diffusion, in Fig.|6]we show 
the theoretical radial distributions in the case when the diffusion co- 
efficient is 10 times less (rather arbitrary) than the value estimated 
for the parameters from Section 3 (we can not put it zero since the 
diffusion term in Equation (2) determines the highest order of the 
derivative). The formation of a gap in the abundance radial distri- 
bution near the corotation in this case is obvious. 

So, as it might have been expected, the main driver for the 
origin of the bimodal radial abundance distribution is the corotation 
resonance. However, in the absence of the turbulent diffusion we 
will have strong depression in radial distribution of heavy elements 
near the corotation. Hence, only combined effect of corotation and 
diffusion leads to a formation of the bimodal distributions both for 
oxygen and iron. The gradient in the inner part derived from iron 
is slightly shallower then the one for oxygen since ^ 1/3 of iron 
is synthesized by long-lived SNe la which are not concentrated in 
spiral arms and are not affected by the corotation resonance. 

In our theory, we aimed to construct a simplest model for ex- 
planation of bimodal radial abundance pattern formation for ele- 
ments, which are synthesized by different sources with different 
galactic distributions. That is why we did not take into account var- 
ious secondary effects like radial inflow of interstellar gas, depen- 
dence of stellar life-time on the metallicity, drift of the corotation 
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Figure 6. The same as in Fig. 3 but for the diffusion coefficient 10 times 
less and = 3.5 kpc. By the arrow is shown the location of the corotation 
resonance. The formation of the deep gap in abundances distributions near 
the corotation in the absence of the turbulent diffusion is obvious. 



resonance, etc. Inclusion of such effects in our model will enable 
to develop a more precise model of formation of the bimodal dis- 
tribution of heavy elements. Also, it should be noted that in the 
current study the proposed mechanism acts together with the ra- 
dially dependent gas infall. However, it will, probably, shape the 
heavy element distribution in a similar manner acting in concord 
with other gradient formation mechanisms, li ke the radially depen - 
dent enr iched galactic wind, a s proposed by IWiebe et al.l ( l200l[) . 
Also, as lAcharova et al.l ( |2005|) have shown, the influence of spi- 
ral arms alone is able to explain the very existence of the gradi- 
ent, not only its bimodal shape. In any case, the discovery of two 
types of SNIa proge nitors (short-lived and long-lived), made by 
iManucci et al. Il l2006h . opens the opportunity to explain not only the 
overall gradient, but the formation of fine structure in radial distri- 
bution over the galactic disc for various heavy elements by means 
of a unified mechanism - the influence of the corotation resonance. 
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